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We formulate a self-consistent field theory for the Chern-Simons fermions to study the dynamical re- 
sponse function of the quantum Hall system at u = ^. Our scheme includes the effect of correlations 
beyond the random-phase approximation (RPA) employed to this date for this system. We report 
results on the density response function, dynamic structure factor, the static structure factor, the 
longitudinal conductivity and the interaction energy of the system. The longitudinal conductivity 
calculated in this scheme shows linear dependence on the wave vector, like the experimental results 
and the RPA, but the absolute values are higher than the experimental results. 
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Despite rapid progress in the field of quantum Hall ef- 
fect in recent years, proper understanding of the state at 
half-filled Landau level still remains a challanging prob- 
lem. A modified Fermi-liquid theory of Chern-Simons 
(CS) fermions, put forward by Halperin, Lee and Read 
(HLR) [Q explained some of the anomalies observed in 
surface acoustic wave experiments (SAW) around the fill- 
ing factor I' = 5 d] • One very interesting result of this 
theory was that at v — ^ the average effective magnetic 
field acting on the fermions vanishes and one expects a 
Fermi surface for those fermions. This result within the 
mean- field approach was later verified in experiments , 
where one finds indications, albeit indirect, for the ex- 
istence of a Fermi surface. In going beyond the mean 
field theory one has to include interactions via the Chern- 
Simons field in order to describe the dynamic response 
functions, transport properties, etc. HLR studied the re- 
sponse functions within the random phase approximation 
(RPA) which takes care of the direct Coulomb interaction 
and the fluctuations in Chern-Simons field. In the work 
of HLR the response functions were analyzed only in the 
long- wavelength limit. The RPA scheme was found to 
explain the wave vector dependence of the longitudinal 
conductivity derived from the SAW experiments ^ . The 
absolute value of the calculated conductivity was, how- 
ever, lower than the experimental results by a factor of 
two. The apparent success of HLR approach has marked 
the beginning of intense activities in the field, so much so 
that often the embellishments tend to overtake the actual 
facts. 

In this letter, we report our studies of the density re- 
sponse function, the dynamic structure factor, the static 
structure factor, and the longitudinal conductivity for the 
quantum Hall system at the filling factor v = where 
we include correlations beyond the RPA scheme for the 
Chern-Simons fermions. In doing that, we have devel- 
oped for the first time a variation on the theme of the cel- 
ebrated self-consistent field theory of Singwi, Tosi, Land 
and Sjolander (STLS) |^ in the quantum Hall regime 
and at a half-filled Landau level. The efficacy of STLS 
approach over the RPA scheme in describing correctly the 



dynamical properties is well established y. Our results 
for longitudinal conductivity show linear wave vector de- 
pendence, as in experiments and also in the RPA scheme, 
but the absolute values are higher than the experimen- 
tal results. Most of the RPA results for the static and 
dynamical functions are also reported here for the first 
time. 

We begin by presenting a few essential steps of the 
HLR approach to establish our notation. The CS trans- 
formation for spinless fermions is defined by (henceforth 
we use units where Ti — e = l) 



§^(r) = *^(r) exp 



-2i / dr' arg(r - r')p(r') 



(1) 



where is the electron creation operator, \E'^'(r) is 

the transformed fermion operator, and arg(r) is the angle 
that vector r forms with the x-axis. The kinetic part of 
the Hamiltonian, which alters due to the transformation, 
is then 



7^kin = :r^ /dr^t(r)[-iV + JA(r)]'*(r) 

2TOfc J 

where is the electron band mass and 
M,(r) = [ dv'cj),{v~v')p{r') 



(2) 



(3) 



[(/)i(r) = 2Vi arg(r)] is the CS field. Expanding the right 
hand side of Eq. (||) and keeping only up to second-order 
contribution, one gets 
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dr^t(r)v2*(r) 



(4) 



X!^^l(^)jkP-k + ^bo(fc) +W2(fc)]PkP- 



where = k x jk is the transverse component of the 
transformed current operator 



j(r)-*^(r) — §(r), 
rrib 



(5) 
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(note that HLR used the diamagnetic current J + 
poA/rrib, with po being the equihbriuni density), where 
Wo = 27r/fc is the Coulomb potential, wi(fc) = 47r/fc, and 
V2{k) = (47r)^/9o/wbfc^ [iejhkhVi{k) is the Fourier trans- 
form of (pj (r)] . This Hamiltonian describes a system with 
the same density as the original system where there is no 
magnetic field but contains a potential vi{k) which cou- 
ples the density fluctuations to the transverse currents. 
This observation is the key ingredient for the exploitation 
of schemes that are normally applied to the electron gas 
in a zero magnetic field. 

We intend to compute the response function matrix 
which gives the density and transverse current responses 
p(k, w) and jXk, w) to external perturbation scalar and 
transverse vector potentials and AS?"^* via 



p(k,cj) 
jT(k,tj) 



where the longitudinal current and vector potential have 
been eliminated using the continuity equation and gauge 
invariance. Following the original derivation of STLS, 
we start from the equation of motion for the one-body 
Wigner distribution function 



(7) 



which determines the density p{r,t) — X]p Z''^'' P! 
and the current j{r,t) = X]p P/''^'' P! the 
semiclassical limit the Heisenberg equation of motion for 
the electron operators at and aj^ gives 

|/W(r,P;i) = ^/(^Hr,P;i) 



rrib 



(V,0,)(r-r')Vp,„ 



-[V,(wo+«2)](r-r')Vp, 



/(2)(r,p;r',p';i) 



+Vp. V,FP"-*(r, t) + Pi Vp,./(i) V.4"*(r, t) (8) 

nib 

where Vp^; = d/dpi, and 
/(2)(r,p;r',p';t) = ^e"'V''''-' 

k,k' 

X («p-k/2 W«P+k/2 W"p'-k'/2 WV+k'/2 W) 

is the two-body distribution function. The key step in 
the STLS approximation consists in the decoupling 

/(2)(r, p; r', p'; t) 2, f('\r, p; t)f^'^ (r', p'; t)g{r - r') 

(10) 

i.e., in the assumption that the correlations in the per- 
turbed, time-dependent state are identical to those in 



the unperturbed, equilibrium state, and are therefore de- 
scribed by the static pair correlation function g(r). No- 
tice that setting g{r) = 1 in Eq. ( p^ ) one recovers the 
RPA where short-range correlations are neglected. 

Equation (^ is equivalent to the response of noninter- 
acting electrons to the effective potentials 



V,FP"-* + / drV(r').9(r-r')V,(«o + i'2)(r-r') 



(11) 



5]jJ(r')5(r-r')V.0,(r-r') 



and 



cff 

j 



Vt j dv'p{Y')g{v-Y')V,4>j{Y-Y') 

(12) 



where the continuity equation has been used to eliminate 
the longitudinal part of the current, and Vt indicates 
projection onto the subspace of transverse vector poten- 
tials. In matrix notation this implies X = x'^il + U'x\, or 
equivalently 



(13) 



where 



X°T 



r ^\ J,' ..0 ) (14) 

is the ideal-gas response function which is known ana- 
lytically |5|[. The matrix of the effective potentials, from 
Eqs. duH^) is 

^^^wo{k) + Mk) ■^Mk)\ 

—iwi(k) ' 



where Wa{k) = [l — Ga{k)]va{k) and the local field factors 
Ga{k) are given by 



G„(fc) = 5][l-5(p)] 



[k.(k-p)] 



(a„+6„)/2 



fc-^" |k-p|' 



(16) 



with oo = 60 = 1; ai = bi = 2, 02 — 0, and 62 = 2, 
and S{k) is the static structure factor, i.e., the Fourier 
transform of the pair correlation function g{r). Using 
the rotational invariance of S{k) it is easy to show that 
in the fc ^ limit Go is linear in fc, G2 is quadratic in 
k, and Gi has a finite limit, Gi(0) = [1 - g{0)]/2 = 1/2. 
Notice that the RPA approximation of HLR amounts to 
Ga{k) = 0, i.e., Wa{k) = Va{k). 

The static structure factor entering Eq. ( p^ ) is ob- 
tained from the fiuctuation-dissipation theorem. 
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Im Xppik,Lu)duj 



(17) 
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FIG. 1. Static density-density response function Xpp(^)O) 
as a function oi k/kp ~ k£o, calculated in the RPA (dashed 
curve) and in STLS (full curve) for — 7. The discontinuity 
in the derivative at = 2fc_F corresponds to the Fermi surface. 



where Xppik,uj)is the density-density response function 
given by Eq. (nl), 



X^pp{k,^) 



1-X°p [wo{k)+W2{k)+wi{kYx''T] 



(18) 



Equations (|13[]1^) are then solved self-consistently for 
a given value of the dimensionless coupling strength, 
Ts — TQ/aB = 2e^/£o^c, where as — l/rribe^ is the Bohr 
radius, ro = (vrpo)^^^^ is the average interparticle spac- 
ing, £o — 1^1^^ is the magnetic length and ojc = B/rub is 
the cyclotron frequency. The relevant values of Vs can be 
estimated in two ways: following HLR ^ we can write, 
rs = 2/C, where C ~ 0.3 is related to the effective mass. 
Alternatively, we can obtain from a realistic estimate 
of ujc and e^/i'o, which is typically, = 1-^3. The 
numerical results show little variation between the two 
cases but for definiteness we consider the first choice. 

Since the density is not affected by the CS transforma- 
tion of Eq. (|l|) , the density-density response function of 
the transformed system is identical to that of the origi- 
nal electron system, and therefore contains information 
about physical properties such as the structure factor, the 
conductivity, etc. In the following we present and discuss 
our results for various quantities derived from Xpp^ ^^'^ 
can therefore drop the distinction between the two sys- 
tems from now on. 

In the static long- wavelength limit (uj ^ kvp, k ^ kp, 
kp = miyVp = l/£o being the Fermi momentmn) one has 



and 



rub 
'2^ 



rrib 



kvp J 



. 2u; 

kvp 



(19) 



(20) 



In the RPA, V2 — pov'j /mb, hence the leading term in the 
denominator of Eq. ( p^ ) cancels, and Xpp(A;,0) vanishes 
linearly for small k. If local field factors are included, 
this cancellation does no longer take place and one gets 
Xpp(fc,0) ~ — (mf,/37r)fc^/fc^. The results are presented 
in Figure |l|, where one clearly sees the difference in the 
limit fc ^ 0. We note that the k^ dependence of the com- 
pressibility at 1/ = 1/2 has been observed also by other 
authors ||^ in the dipole nature of = 1/2 state, which 
arises primarily due to projection to the lowest Landau 
level. Further, we find that excluding the cyclotron con- 
tribution, /Im Xppwdw oc and /Im Xppdw oc apart 
from possible logarithmic terms. 

The longitudinal conductivity, which is relevant to 
surface-acoustic-wave experiments is given by <j~]: — 
{{k"^ /uj)[x~l^p{k,(ju) — Xpp^(A:,0)]. Since the speed of sound 
Cs is small compared to the Fermi velocity vp and k <C 
kp, we can use the limiting forms (p^po|) into Eq. ([T^). 
This leads to the result a^xik, Cgk) ~ k/2Trkp, which 
has the same linear dependence on fc, but is twice the 
experimental values The RPA result of HLR is 

'^xx^{k,Csk) ~ k/Snkp. A quantitative agreement with 
experiment can however be achieved if the CS interaction 
4>j{r) is softened at small separation. 

Using the well-known asymptotic behaviors, x'^pp — 
pok'^/nibUi^ and Xt — 0{k'^), valid for k£o <^ 1, kvp <^ uj, 
one sees that the dynamical structure factor 



S{k,uj) = —IiLnxpp{k,uj) 



(21) 



has a pole at the cyclotron frequency u; = lUc, describing 
inter-Landau- level excitations, which - at fc = - is un- 
affected by correlations, and is in agreement with Kohn's 
theorem. The mode dispersion, which is computed by 
locating the zeroes of the denominator in Eq. (|8|) , turns 
out to be significantly lower than the RPA result (see 
inset of Fig. Our finite-fc results are presented in 

Figure ^ where we have a (5- function peak at w ~ Wc, 
corresponding to the cyclotron motion and a continuum 
of particle-hole excitations in the range 



i.2 i.2 
vpk < \uj\ < 



2mb 



2mb 



■ Vpk . 



(22) 



Our results for the static structure factor S{k) are plot- 
ted in Fig. 1^. Here we compare our RPA results, calcu- 
lated from Eqs. (|l^|l|) with Wa = Va, and the STLS 
results, at = 7. These results are also compared with 
S{k) calculated for a modified Laughlin state a.t v = ^ 
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FIG. 2. Dynamic structure factor S{k,w) for k = O.Gkp, 
as a function of lj (in Rydberg) in RPA (dashed curve) and 
STLS (full curve). The 5-function peak corresponding to the 
inter-Landau-level mode has been artificially broadened for 
clarity and contains most of the spectral strength. The inset 
shows the excitation spectrum, composed of the particle-hole 
continuum plus the sharp cyclotron mode. 



, proposed by Read |g] . All these curves obey the lead- 
ing [kiof /2 behavior at small k. As expected, the STLS 
scheme includes substantial amount of correlations and 
hence is significantly higher than the RPA results near 
k = 2kF- 

Knowledge of the structure factor S{k) allows us to 
compute the potential energy per particle, {PE) — 
\ VQ{k)[S{k) — l]. The full interaction energy (defined 
as the total energy minus the noninteracting term hujc/'^) 
is then obtained via coupling constant integration. 



1 



dr's{PE{r',)) 



(23) 



(we measure energies in units of e^/^o). The STLS result 
i^int — —0.48 compares favorably with finite-size exact di- 
agonalization studies i^int = —0.466. The RPA over- 
estimates appreciably the interaction energy, and gives 
iJint — —0.76. At the same rs, the STLS potential energy 
is {PE) = —0.49, showing that the inter-Landau-level ki- 
netic energy is a minor contribution. 

In summary, we have presented a self-consistent 
scheme for the calculation of the dynamical response 
function of a quantum Hall fluid at i/ = i, based on a 
generalization of the STLS method to the case of Chern- 
Simons fermions. Our results exhibit significant differ- 
ences with the RPA computations, in particular on the 
longitudinal conductivity, the static response function 
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are also compared with the results from Ref. 



3. Static structure factor S{k) as a function of 
kio in the RPA and the STLS scheme. The results 

" (TC). 



and the structure factor. 
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